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In the strictly periodic setting, the electric polarization of inversion-symmetric solids with and 
without time-reversal symmetry and the isotropic magneto-electric response function of time-reversal 
symmetric insulators are known to be topological invariants displaying an exact Z2 quantization. 
This quantization is stabilized by the symmetries. In the present work, we investigate the fate of such 
symmetry-stabilized topological invariants in the presence of a disorder which breaks the symmetries 
but restores them on average. Using a rigorous analysis, we conclude that the strict quantization still 
holds in these conditions. Numerical calculations confirm this prediction. 

PACS numbers: 73.43.Cd,73.43.Nq,74.62.En,75.85.+t 


An important issue in the field of symmefry-profecfed 
topological phases is fhe fafe of fhe topological invari- 
anfs in fhe presence of bulk and surface disorder which 
can break fhe symmefries 111. For fhe fime-reversal sym- 
mefric topological insulafors, fhis issue is relafed to fhe 
confaminafion of fhe samples wifh magnetic impurities 
and if has been addressed experimenfally and fheoref- 
ically in a number of works m. While fhis confami¬ 
nafion can be fully confrolled in laborafory condifions, 
fhis is nof fhe case in fhe real world condifions. Still, 
if fhe confaminafion is small, fhe magnetic impurities 
are in a non-correlafed disordered phase and, on av¬ 
erage, fhe fime-reversal symmefry is preserved. Then 
a question which is is exfremely relevanf for fhe prac¬ 
tical applications of fhese materials is if fhe fopologi- 
cal characferisfics, such as fhe exfended characfer of fhe 
surface sfafes and fhe quanfizafion of fhe bulk fopologi- 
cal invarianf, are preserved under such "average" fime- 
reversal symmefry condifions? The firsf characferisfics 
has been shown in Ref. ||3l fo hold if fhe disorder is nof 
too sfrong and here we show fhaf fhe isofropic magnefo- 
elecfric response function 21 remains quantized in such 
disordered regimes, provided a specfral gap is presenf. 
The anfiferromagnefic fopological insulafors infroduced 
in Ref. |[5| is anofher class where fhe topology is sfa- 
bilized by an average fime-reversal symmefry. In fhe 
presence of disorder, fhis class of insulafors were re- 
cenfly shown ||6l fo posses disfincf fopological phases 
surrounded by sharp phase boundary which can be de- 
fecfed by fransporf measuremenfs. Furfhermore, fhe 
weak fopological insulafors can be fhoughf as profecfed 
by fhe franslafional symmefry and disorder definifely 
breaks fhis symmefry buf fhe symmefry is resfored on 
average. Refs. IZH^ gave evidence fhaf, in cerfain con¬ 
difions, fhis is enough for fhe surface sfafes of weak 
fopological insulafors fo avoid Anderson localization. 

When some of fhe sfabilizing symmefries are space 
symmefries, such as for fhe crysfalline fopological in¬ 
sulafors in |T0l, fhe inversion-symmefric lITTl [T2ll . fhe 
reflecfion-symmefric ITSl [141 or fhe spin-orbif free fopo¬ 
logical insulafors IT^ in general, fhese issues are acfually 


of cenfral imporfance because fhe space disorder of fhe 
afoms is inevifable and impossible fo confrol even in 
fhe laborafory condifions. If, for example, one consid¬ 
ers fhe random displacemenfs of fhe afoms due fo finife 
femperafure, fhere can be no expecfafions fhaf fhe dis¬ 
ordered pofenfial respecfs fhe underlying symmefry of 
fhe lattice buf one can be sure fhaf, on average, fhe sym¬ 
mefries are sfricfly preserved since fhe fhermod 5 mamic 
sfafe (including fhe nuclei) is symmefric. Such resfora- 
fion of fhe symmefries by averaging also occurs when 
fhe disorder is induced by defecfs which do nof desfroy 
fhe crysfalline order complefely. Ref. Ill infroduced 
a Z 2 fopological invarianf, which in principle covers 
all classes of symmefry-profecfed fopological insulafors 
(and more) in such disordered condifions. This fopolog¬ 
ical invarianf was puf fo work for a disordered model 
wifh averaged reflecfion symmefry and a localizafion- 
delocalizafion fransifion was observed numerically ex- 
acfly af fhe poinf where fhe invarianf changed ifs value. 
In fhe presenf work, we consider disordered insulafors 
wifh averaged inversion symmefry and show fhaf fhe 
classical elecfric polarization assumes sfricf quantized 
values 0 or provided a specfral gap is presenf. 

Now, fhe exisfing definifions of fhe symmefry- 
sfabilized fopological invarianfs depend fundamenfally 
on fhe exactness of fhe symmefries, hence, af fhe firsf 
sighf, if seems impossible fo define fhem for a concrefe 
disorder configuration which, jusf by ifself, breaks fhe 
symmefries. However, as noted in Ref. m, when fhe 
symmefries, which are preserved only average, are com¬ 
bined wifh fhe franslafions fhen a cerfain self-averaging 
properfy fakes place, enabling one fo define exact fopo¬ 
logical invarianfs for such realistic condifions. In fhe 
presenf work, we provide a precise formulation of fhis 
self-averaging properfy wifhin fhe framework of homo¬ 
geneous disordered sysfems ll^riTl . While af fhis mo- 
menf we carmof make any sfafemenfs abouf fhe bound¬ 
ary sfafes, our findings definifely confribufes fo fhe 
growing body evidence fhaf fhe symmefry-sfabilized 
fopological insulafors are more robusf fhan previously 
fhoughf. 
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The paper is organized as follows. In the first sec¬ 
tion, we describe the homogeneous disordered systems 
and we formulate a precise relation between the disor¬ 
der and symmetry, called here the compatibility con¬ 
dition, which ensures the restoration of the symmetry 
on average. An explicit yet very generic model for 
homogeneous disordered solids is presented, and the 
inversion as well as the time-reversal symmetries are 
shown to be compatible with the model. The second 
section discusses the electric polarization of homoge¬ 
neous disordered systems with inversion symmetry and 
no time-reversal symmetry. The disorder is assumed 
to be compatible with the symmetry. The polarization 
is shown to be an intensive macroscopic function with 
a self-averaging property (preventing fluctuations from 
one disorder configuration to another), and which takes 
only quantized values. The fourth section reports a nu¬ 
merical analysis which confirms the theoretical predic¬ 
tions. The fifth section applies the arguments to the 
isotropic magneto-electric response function of TRS in¬ 
sulators, leading to similar conclusions. Last section 
summarizes our conclusions. 

HOMOGENEOUS SOLIDS 

Physically, a homogeneous system is defined as an ex¬ 
tended system with translation invariance broken at the 
microscopic scale but this symmetry-breaking is unde¬ 
tectable at macroscopic scales. In other words, taking 
micron-size samples from a large piece of a homoge¬ 
neous material will lead to same intensive TD parame¬ 
ters and coefficients (of course not exactly the same for 
finite pieces, but the differences are well below the ex¬ 
perimental resolution). Formulating this property in a 
mathematically rigorous way was quite a challenge for 
the mathematical physics community but these days the 
concept of a homogeneous solid state system has a very 
precise meaning and the mathematical framework built 
around this concept is natural and fruitful for a large 
class of aperiodic systems IfThlllTl . According to the pre¬ 
cise definition, an aperiodic tight-binding Hamiltonian 
H over the lattice is homogeneous if H together with 
its lattice translates TaHT~^, for all a e define a set 
which has a compact closure in the strong topology of 
bounded operators. Here and throughout, will de¬ 
note the lattice translation by a. What we are going to 
present in the following applies strictly to the homoge¬ 
neous system so defined, but for concreteness we will 
carry the discussion in the context of disordered lattice 
models, which are explicitly constructed in Example 1. 

Henceforth, let us consider a generic finite-range dis¬ 
ordered lattice model 

H^=Ho + ( 1 ) 

defined over the Hilbert space spanned by \n, a), where 


n e represent the nodes of the lattice indexing the 
unit cells and a = 1,..., Q the orbitals associated with 
a unit cell. The orbital index includes the spin de¬ 
gree of freedom, but the latter will be separated out 
when needed. The orbitals can and will be chosen 
to be real, that is, invariant to complex conjugation: 
'K\n,a) = \n,a). The translational invariant piece Hq 
of the Hamiltonian is assumed to depend on a set of 
N parameters 4 = (^i/ • • ■ / 5 n)- When needed, we write 
this dependence explicitly as Hq{^) or H^{E,). We include 
such dependence because both the electric polarization 
and the magneto-electric response functions are defined 
through deformations of the system from a reference 
configuration. If the reader is more comfortable with an 
explicit representation of this parameter space, he may 
think of the vector ^ as the (always finite) collection of 
hopping amplitudes. 

The random potential depends on the disorder 
configuration co, which is seen here as a point in a 
disorder configuration space Q, which is compact and 
equipped with a probability measure dco. The system 
is assumed to be homogeneous, which in this context is 
assured by the covariant property 

TaH„T-^ = Vfl e Z^ (2) 

where t represents a homeomorphic action of the group 
of lattice translations on Q. The action t is assumed 
ergodic and probability preserving so that Birkhoff's er- 
godic theorem IITHI applies, namely, the following iden¬ 
tity holds with probability one in o): 

[ dco'ficv'). (3) 

It is precisely this identity which ensures the non¬ 
fluctuating character of the intensive thermod 5 mamic 
functions from one disorder configuration to another. 

We now formulate a precise condition which automat¬ 
ically leads to the restoration of the symmetry after the 
disorder average is taken. Consider a generic symmetry 
operation: 

S\n,a) = y^R„a-,„'a'\n ,a'), (4) 

n',a' 

which can be linear or anti-linear. Recall that the orbitals 
are real, hence the symmetry is fully determined by the 
coefficients Rna;n’a’- The symmetry generates an action 
on the parameter space: 

(5) 

and on the disorder configuration space: 

( 6 ) 

for which we use the same notation S. As a consequence, 
we have the following covariant property: 


(7) 
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In fact, such covariant property is obeyed by any func¬ 
tion of the Hamiltonian: 

S0(h^{^))S-^ = (8) 

Definition. We say that the symmetry is compatible with 
the disorder if 

d{Sco) - dco, (9) 

that is, if the probability measure dco is invariant under 
the action on Q induced by the symmetry. 

Perhaps the reader already noticed that we are merely 
restating the condition formulated in Ref. m Section IIA, 
in a more concise mathematical notation. This definition 
is relevant when we consider the disorder average: 

5 dco o(h„(^))J5-1 = f^dco 0 {Hso,(S^)), (10) 


which leaves dco invariant and is known to act ergodi- 
cally, hence Eq. applies. It is straightforwad to check 
that the Hamiltonian has indeed the covariant property: 

TaHa^T-,^ = (17) 

Example 2: Compatibility of the inversion symmetry 
with the model^^ This symmetry maps the unit cell n 
into -n and it can mix the orbitals in the process. How¬ 
ever, we can always choose an orbital basis so that no 
such mixing occurs, and since - 1: 

I\n,a) = Y^Oa\-n,a), (18) 

a 

where all a's are signs. The induced action on Q is: 

{Icofnt = OaOliC0‘%^-,„. (19) 

Now, 


because a change of variable Sco co leads to: 


r d{S-^co) o(h„(5^)) = f dco 0 (hUS^). (11) 


dco^ Yl = n n ‘^'^tmdcO^Zn (20) 

n,m)a n,m;a<^ 

and now one can see explicitly that d{Ico) = dco. 


In other words 

5 (^£ dco o(h,,(^))| 5-1 = £ dco o(h,,(5^)). (12) 

In particular, if ^ is a fixed point for the symmetry, 5^ = 
4, or equivalently if SHo(4)S~^ = Hq{4), then 

5 dco o(h,,(^)J5-1 - £‘ico ^(hU4)\ (13) 

which shows that Eq. implies the restoration of the 
symmetry for the averaged quantities. 

Example 1: Generic homogeneous disordered model. 

£ {C_J4) + Wco:Q\n,a}(m,fil (14) 

where are independent random entries drawn uni¬ 
formly from the interval [-|, |]. The collection of all 
random variables co = {co‘^^,„} can be viewed as a point in 
an infinite dimensional configuration space Q, which is 
just an infinite product of intervals [-|, |]. The result is 
a compact and metrizable Tychonov space which can be 
equipped with the product probability measure: 


Example 3: Compatibility of the time-reversal symmetry 
with the model\l^ The TRS is defined as the anti-linear 
map 


0|M,a,cT) =-cr|M, a,-u), (21) 


where the spin degree of freedom cr = +1 (for spin 
up/down) was separated out. Note that 0^ = -1. The 
induced action on Q is: 


Then, 


(0(U) 


a,0)^,cy' 

nm 


oa'co 


n,m 


( 22 ) 


dco = jQ 


dco. 


a,o;^,o' 


(23) 




n 


a,l;;5,l -1 


j j I 

dco„C dco. 


1 J u 

aco„ 


and now one can see explicitly that d{&co) = dco. 


ELECTRIC POLARIZATION OE INVERSION 
SYMMETRIC INSULATORS 

Generic definition 


dco= Yl 

n,a;tn,^ 


(15) 


There is a natural action of the lattice translations on Q: 


{iaCO)„% = (On-a,m-a> 


By definition IIT9ll20l , the change in the electric polar¬ 
ization, as a result of a macroscopic deformation Ha,(4t) 
of the Hamiltonian, is: 

= f ja,{f) dt, 

Jo 


a e lP, 


(24) 
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where ja,{t) is the density of the charge-current. The 
latter is microscopically defined as: 

j Q 

ja,it) = lim - y y\(n, a\pco(t)Ja,(t)\n, a), (25) 

V—»oo V * ■ * ^ 
neV a=l 

where pa){t) is the time-evolved density matrix and 

Ut)-ie[HU^t),X] (26) 

is the current operator. Here, X denotes the position op¬ 
erator and, for convenience, fhe elecfron charge e will 
be sef fo unify in fhe following. If fhe deformafion 
sfarfs from fhe fhermodynamic equilibrium sfafe, fhen 
fhe fime-evolved densify mafrix is 

put) = Ut(^FD{HU^o))u;\ (27) 

wifh Ut being fhe unifary time evolution generated 
by HU^t) and Ofd the Fermi-Dirac distribution. One 
should not confuse pUt) and <l>FD(hfa,(^t))/ because fhe 
fime-evolved densify mafrix is no longer given by fhe 
Gibbs sfafe. 


Self-averaging property 


Now nofe fhaf fhe fundamenfal formula in Eq. [^in¬ 
volves fhe densify of fhe currenf rafher fhan fhe currenf 
ifself, and fhis is why fhe frace per volume appears in 
Eq. 25 This is consisfenf wifh fhe facf fhaf polarizafion 
is an intensive macroscopic function. Erom a technical 
point of view, fhis is an imporfant observation because 
we can use Birkhoff's ergodic theorem ITHl to demon¬ 
strate that ISfPa, is, with probability one, independent of 
the disorder configuration co. Indeed, note that the op¬ 
erator inside the brackets of Eq. 25 called in fhe lines 
below, is covarianf: 


TaFo,T-^ = (28) 

Then: 

i Y,<n\FUn) = F Y^{0\F,-,J0) = f dco (0|F,„|0) (29) 

neV neV ^ 

in fhe limif V —> oo. The conclusion is fhaf AP^) is self¬ 
averaging and ifs macroscopic value comes as an average 
over fhe disorder configuration space. As such, we can 
drop fhe subscripf tu in APa> from now on. 


Schulz-Baldes-Teufel formula 

We consider now infinifely slow deformations of fhe 
Hamiltonian, which are better visualized as pafhs y in 


fhe paramefer space, paramefrized as !^s}s€[o,i]- By em¬ 
ploying fhe adiabafic fheorem, Schulz-Baldes and Teufel 
showed in Ref. IZDl fhaf in fhe exfreme adiabafic limif 
and when fhe femperafure goes fo zero: 

A!P(y) = dsT(PU^s)i[dsPU^s)A[X,PU^,)]]), (30) 

where 7~ denotes fhe frace per volume: 

Q 

7'(...) = lim y<n,a|...|n,a), (31) 

V—»oo / J / J 
n€V a=\ 


and 

PUE,) = y[-oo,.q(H,,(^)) (32) 

is fhe Fermi projection onto fhe occupied sfafes af coor¬ 
dinate . 5 . Above, if is assumed fhaf fhe specfral gap of 
remains open for all s e [0,1] and fhaf fhe Fermi 
level £f is inside fhis gap. We will keep fhe i = 
in fronf of fhe commufafors in order fo make fhem self- 
adjoinf. Now, using Birkhoff's fheorem as before, we can 
wrife, equivalently: 


AP{y) - f ds f dco 

Jo Ja 

{0,a\PU^s)i[dsPU^s), i[X, PU^J]]\o, a). 

This is Schulz-Baldes-Teufel formula, which can be re¬ 
garded as fhe disordered version of fhe King-Smifh- 
Vanderbilf formula for fhe sfafic sponfaneous orbifal po¬ 
larizafion 1221 . 


Q 

E 


Quantization 


The proof of quanfizafion proceeds in fwo sfeps. Firsf, 
based on Schulz-Baldes-Teufel formula Eq. we can 
demonsfrafe fhaf fhe change in fhe elecfric polarizafion 
along the inverted path Jy (see Fig.j^ is: 

AP{Iy) = -AP{y). (34) 


This equality is remarkable because we are computing 
the polarization of a system in an arbitrary disorder con¬ 
figuration which breaks the inversion symmetry. For the 
periodic case with strict inversion symmetry, this prop¬ 
erty is well known IfTTl [T^ . The keys to its proof are 
fhe self-averaging properfy fhe compafibilify befween 
fhe inversion symmefry and disorder, and fhe behavior 
of fhe position operafor under inversion, IXI~^ = -X. 
The proof proceeds as follows: 

AP{Iy) = r ds f dco 
Jo Ja 

(O, alPUI4s)ilJsPUJ'U i[X, F„( J4)]] [O, a). 


Q 

E ™ 
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and observe that 


(36) 


Then a change of variable co Ico and cancelations of 
ferms like II~^, lead us fo: 


J r^l ^ Q 

ds d{Ico) y 

0 


(37) 


(O, a\lPU^s)i[dsPa,{^s), i[-X,P„(4,)]]-r“^l0, a). 

The probabilify measure is invariant, d{Ico) = dco. Also, 
the 0-site is left invariant by I and, since we are tracing 
over the orbital degrees of freedom, fhe acfion of fhe re¬ 
maining I operators have no effecf and can be removed. 
Then Eq.[34|follows. 


The usual consfanf in fronf is absenf here because one 
derivafive is faken in fhe 'k'-space and fhe ofher in fhe 
real-space represenfafions. As tongs as Pq,{^) is smoofh 
along fhe loop and ifs mafrix elemenfs {n\Pa,{^\m) decay 
sufficiently fast, all these Chern numbers are integers for 
dimension D = lor2. InD = 3 fhese Chem numbers 
are weak topological invariants and they remain integer 
only if magnetic fields are nof presenf. All fhese condi- 
fions are mef if fhe specfral gap af fhe Fermi level remains 
open along fhe loop. 

In conclusion, when insulafors are deformed befween 
fixed poinfs of fhe inversion symmefry operafor, fhe 
change in fhe elecfric polarization is quantized in units 
of even fhoughf fhe symmefry exisfs only for disorder 
averages. 



FIG. 1. The direct path y (upper semi-plane) and the inverted 
path -ly (lower semi-plane) in the parameter space The 
horizontal axis represents the manifold which left invariant by 
the symmetry. The interesting cases occur when the spectral 
gap closes along this manifold, which is schematically shown 
by the gray region. 


The second step of fhe proof consisfs of fhe following 
argument. Assume that the initial and final Hamilfoni- 
ans are inversion symmefric, fhat is, and are fixed 
poinfs of I, 1^0 - E,o and = E,i. Then any pafh y 
joining fhe two Hamiltonians can be closed info a loop 
by augmenting -ly, where fhe minus sign imply fhaf 
fhe pafh is walked in reverse (see Fig. 1). Nofe fhaf fhe 
argumenf doesn'f work if fhe end poinfs of y are nof 
fixed poinfs for I. Given Fq. 34 fhis has fhe following 
effecf: 

h!P{y) = i/ 2 (A!P(y) -h ^'P{-Iy)) = ^liMPiy - ly). (38) 


Numerical confirmation 


We consider the 1-dimensional Rice-Mele model 
«-(^) =5 E + (-l)".Si)(|nXn + 1| + In + l)(n|) (41) 
+ \Y,((-\)%2 + V^ayn) \n){n\, 


with onsite disorder. As detailed in 11261 , this model was 
introduced originally by Rice and Mele to study solitons 
in conducting polymers but it was subsequently used in 
other contexts too. It is a particular form of fhe generic 
model Fq. 14 hence everyfhing we sfafed abouf Fq. [T4| 
continue fo hold here. In fhe absence of disorder, fhe 
Rice-Mele model is gapped if half-filling is assumed, ex- 
cepf af ^ = 0, and if is inversion-symmefric whenever 
^2 = 0. Also, if is well known 12^ fhaf adiabatic de¬ 
formations around closed loops surrounding fhe origin 
leads fo a quanfized polarization change hfP = +1 (or 
Cl = +1). 

The model has fwo sfafes per unif cell and if can be 
brought to the form wriffen in Fq.j^by defining 


|2n) = |n,l), |2n-f-1) = |n,-l). 


(42) 


Buf as already noted in 11211123] , fhe change in fhe polar- 
izafion along a closed loop leads fo fhe non-commufafive 
Isf Chern number defined in Ref. 124] . The conclusion 
is: 


APXt) = V2 Ci[(y - ly) X Sj], ; = 1,..., D, (39) 

where on fhe righfhand side is fhe Chem number of fhe 
Fermi projecfion over fhe manifold (y — ly) times fhe 
secfion of fhe non-commufafive Brillouin forus along fhe 
jih direction: 

Ci= f d^T(P^(^)i[diP^(4).ilXj,P^(^)]]). (40) 

y-Jy 


in which case fhe Hamiltonian fakes fhe form: 


n€Z a=±l 

X (jn,aj(n + - a),-a\ + \n + f(l - a),-aj(n,a\'j 

HEE (a 52 + fVcu„,a;) |n,a)(n,a|. 

M€Z a=±l 


If in fhe original rendering of Fq. 41 we choose fhe 
inversion point to be between the 0-th and first site, then 
the inversion operation takes form: 


J|n,a) \ -n, -a). 


(44) 







6 






(a) 


CL 

CO 

CD 

"cO 

o 

CD 

CL 

<X> 


0.8 

0.6 

0.4 

0.2 

O 


111 

-1-1-1-1- 

(a’) 

-1 





- 1, 





\ 

iilllliiilili 

if 

mil II 

lii'iiii 

1 1 1 1 1 1 1 



1234567 

Disorder strength 



1 

1 1-1 1 1 


-■ 


III 1 

iiiiiiiii 


.mill 

. . 










. i ; 



(C) 

i j_i i 1_1 



Disorder strength 


Disorder strength 


FIG. 2. (Color online) Top row: The variation of the electric polarization when the Rice-Mele model is adiabatically deformed 
along a semi-circle (black dots) and along a full circle (red dots) in the two-parameter space. In the latter case, the variation is 
equal to a Chern number, hence the label. The calculations are performed in the presence of a disordered potential which breaks 
the inversion-symmetry but restores it on the average. The disorder strength was increased from 0 to 7. Each dot corresponds 
to one disorder configuration and 100 disorder configurations were considered for each disordered strength. The average and 
statistical variance are not shown on purpose (see text). Bottom rows: The spectral gap of fhe Hamiltonian at half-filling. A 
dot represents the spectral gap at one disorder configuration and one position along the adiabatic deformation path. The three 
columns correspond to the different system sizes: L = 100 (a), 200 (b) and 400 (c). 


It can be set in the diagonal form of Eq.j^if we choose 
to work with the states 1) ± n, -1)). However, we 

will not do that here. A direct computation will show 
that: 

= Hra,(IE,), (45) 

with 

^ 2 ) = (^ 1 , -«2), {Ia))n,a = CO-n.-a- (46) 

Let us stress again that the onsite disorder breaks the 
inversion symmetry because we impose no correlation 
between co„^a and co-„-a- However, the probability mea¬ 
sure Yin,a dc0n,a is easily seen to be invariant under I. 
We now choose the path y to be the semi-circle: 


y = 0.1 X (coss,sins}se[o,7i]. 


(47) 


which connects two fixed points of the inversion sym¬ 
metry. By augmenting with -lywe obtain the full circle, 
s e [0, 27 i]. The numerical calculations are performed af 
finife size n e {0,..., L} and wifh periodic boundary con¬ 
ditions. Formulas from Eqs. 30 and 40 are used for all the 
results reported here. The commutator with the position 
operator is implemented using the strategy developed 
in Ref. |[27l . The pafh infegral is discretized using 100 
poinfs for y and 200 poinfs for y — ly. The derivative 
wifh respecf fo 4 along fhe pafh is compufed using fhe 
five-poinf sfencil fmife-difference approximafion. The 


disorder sfrengfh W was increased gradually and fhe 
calculation was repealed 100 times wifh updafed ran¬ 
dom pofenfial for each W. 

The resulfs are reporfed in Fig. 2, where each dof seen 
in these plots represents a single disorder configuration. 
In other words, no disorder average has been performed on 
the data. In the first row we show the results for A!P(y) 
and Chem number or A!P(y - Iy). The firsf fhing fo 
notice is fhe clusfering of fhe dafa for A!P(y) around fhe 
predicfed quanfized value of af disorder sfrengfhs 
lower fhan a critical value W^.. ITae flucfuafions around 
fhe quanfized values can be affribufed fo fhe finife sys- 
fem sizes because fhe flucfuafions can be seen fo dimin¬ 
ish as fhe size of fhe sysfem gefs larger. This confirms our 
fheorefical predicfion fhaf A!P(y) is self-averaging hence 
non-flucfuafing in the thermod 5 mamic limit, and that, at 
least for small W's, if fakes quanfized values in fhe unif 
of 

The Chern number shows a very good quantization, 
wifh virfually no flucfuafions from one disorder configu¬ 
ration to another below the same W^. The reason for fhis 
difference is fhaf fhe quanfizafion of fhe Chem number 
does nof require fhe restoring of fhe symmefry. If only 
requires fhat fhe exponenfial decay rate of fhe Fermi pro- 
jecfion be large compared fo 1/L. The Fermi projecfion 
also needs fo be smoofh of E, along fhe loop, buf fhis is 
aufomafically fhe case if fhe specfral gap remains open 
(which is fhe case for small W's). If, for example, we av- 
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erage A!P(y) over just two disorder configurations which 
are mirrored by inversion symmetry, then the quantiza¬ 
tion of A'P{y) will be as good as fhaf of fhe Chem num¬ 
ber (fhis can be shown exacfly). While fhis could be a 
beffer numerical mefhod fo compufe A'P{y), if is irrele- 
vanf for fhe presenf discussion because here we wanf fo 
demonsfrafe a principle, namely, fhaf fhe quanfizafion 
of A'P{y) occurs even for a single disorder configurafion 
CO, provided fhe sysfem size is large enough. So, whaf 
is acfually happening in fhis laffer case? For fhis, lef 
us nofe fhaf fhe formula for A'P{y) remains invarianf if 
fhe disorder is franslafed. Now, due fo fhe ergodicify of 
fhe franslafions, when we franslafe fhe disorder we sfarf 
exploring fhe disorder configurafion space Q, and if fhe 
size of fhe sysfem is large enough, af some poinf a frans- 
lafion will bring us close fo a new disorder configurafion 
which looks in many respecfs like fhe inverfed disorder 
1(0. In ofher words, if fhe sysfem size is large enough, 
we don'f need fo average over two mirrored disorder 
configurafions because fhis already happens due fo fhe 
ergodicify of fhe sysfem. 

The lasf issue we need fo address is fhe value of fhe 
crifical W^. We know already fhaf Wc should be larger or 
equal fhan fhe value of W where fhe specfral gap of (^) 
closes somewhere along fhe loop. Can W^ be sfricfly 
larger fhan fhis value fhaf we menfioned? To answer 
fhis quesfion, we ploffed in fhe second row of Fig. 2 fhe 
specfral gaps of af half-filling, as ^ was varied 

along y and fhe disorder configurafions were updafed. 
These many dafa were fhen collapsed and shown for 
increasing values of W. Whenever a dof in fhe second 
row of Fig. 2 touches fhe horizonfal axis, fhe specfral 
gap of closed for some 4 along fhe loop and some 

disorder configurafion. Nofe fhaf fhe flucfuafions of fhe 
gap due fo disorder die ouf in fhe fhermod 5 mamic limif. 
By comparing fhe firsf and second rows of Fig. 2, we can 
conclude wifh confidence fhaf fhe deviafions from fhe 
quanfized values in fhe fop row of Fig. 2 occur exacfly 
af fhe value of W where fhe specfral gap closes. The 
answer fo our quesfion is no, and fhe reason is because 
fhe Fermi projecfion fails fo be smoofh of ^ beyond fhaf 
poinf. 


MAGNETO-ELECTRIC RESPONSE OE TRS 
INSULATORS 


fhe argumenfs are repefifive, we will expedife fhe expo- 
sifion. 

The main fool of our analysis is fhe formula derived 
in Ref. 1281 for fhe change Aa during a deformation of 
fhe sysfem along a pafh y = {^s}s€[o,i] in the parameter 
space: 

Aa{y) = dsT(Pa,{Qd,,PUQ ■ ■-d^^PU^s)) 

where 7* denotes the trace per volume, dj is a shorthand 
for: 


djPU^s) = i[Xj, PUQl for 7 = 1,2,3, (49) 


and diP,oi^s) = dsPioi^s). Also, is fhe anfisymmef- 
ric fensor and summation over repeating indices is as¬ 
sumed. Since fhe operafor inside fhe frace per volume is 
covarianf, we can apply Birkhoff's ergodic fheorem IITSlI 
fo equivalenfly write 


Aa{y) 



(50) 


y, (o,a|P,y(^s)^!iF’„(4). ..di^P,„{^s)\o, ay 


This shows af once fhe self-averaging properfy of fhe 
magnefo-elecfric response funcfion. 

Nexf, we show fhaf if fhe deformation occurs along 
fhe TRS mirrored pafh 0y, fhen: 


Aa(©y) = -Aa(y). (51) 

This properfy is well known for periodic ID and dis¬ 
ordered Il28l TRS insulafors, buf here we compufe fhe 
magnefo-elecfric effecf for a sysfem in a disorder con¬ 
figurafion which breaks fhe TRS. The proof proceed fhe 
same way as for polarization, using fhe self-averaging 
properfy, fhe compafibilify between TRS and disorder, 
together with the behavior of the derivatives under TRS: 
0(9y0“^ = —dj, for j = 1,2,3, and &ds&~^ = ds. The proof 
goes as follows: 


Aa(0y) =le'i •■‘4 



(52) 


y {o,a\PU@Qd„PU®^s) ■ ..d,,PU®Q\0,a). 


Consider now an insulator H,y(^) in dimension D = 3. 
In this section we investigate the isotropic part of the 
magneto-electric response function: 



dPj 

~Wj' 


(48) 


If TRS is considered, then the partial derivatives with 
respect to the magnetic field are taken at B = 0. Since 


Recall that 


Pa,(@^s) = @Pe-UQ®~^. (53) 

Then, after a change of variable co —> &co: 

... = - f ds f d{&co)y (54) 

Jo Ja a 

(O, ajePUm^PMs) ■ ■ ■ dM^s)@-yo, a). 
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Then the statement follows from fhe explicif acfion of 0 
given in Eq. [^and fhe compafibilify befween TRS and 
disorder, d{&co) = dco. 

Lasfly, given a pafh y befween fwo TRS fixed poinfs, 

= 0^0 and 4i = 0^i, one can close fhis pafh info a loop 
by augmenfing wifh ifs TRS mirrored image -@y, faken 
wifh opposife orienfafion. Nofe fhaf fhis argumenf will 
nof work if fhe end-poinfs of y are nof TRS fixed poinfs. 
Then: 

Aa{y) = yi^Aaiy) + Aa(-0y)j = ^/ 2 Aa{y - 0y). (55) 

For fhe periodic TRS case, if is well known IHl fhaf 
fhe variafion of fhe magnefo-elecfric response funcfion 
along a closed loop is equal fo a second Chem number. 
As demonsfrafed in Ref ll28l , fhis remains frue in fhe dis¬ 
ordered case, in which case fhe connecfion is wifh fhe 
non-commufafive second Chem number infroduced in 
Ref. 12^ . More precisely: 

Aa(y - 0y) = CzfCy - 0y) x T^], (56) 

where on fhe righf we have fhe second Chem number 
of fhe Fermi projecfion Pai{^) over fhe manifold (y - 0y) 
fimes fhe 3-dimensional non-commufafive torus: 

C2 = e'-‘* f (57) 

y-0y 

This is a sfrong fopological invarianf which is known 
Il29l fo fake only infeger values. Nofe fhaf fhe consfanf 
in fronf differs from fhe usual consfanf because some of 
fhe derivafives are faken in k-space and some in real- 
space. 

In conclusion, when 3-dimensional insulafors are de¬ 
formed befween fixed poinfs of fhe fime-reversal oper- 
afion, fhe change in fhe isofropic magnefo-elecfric re¬ 
sponse funcfion is quanfized in unifs of j, even fhoughf 
TRS applies only for disorder averages. Recall fhaf fhe 
exisfence of a specfral is required by our argumenf. The 
numerical simulafions of Aa performed for fhe work 
Ref. Il30l indicafe fhaf, as in fhe previous case, fhe quanfi- 
zafion does nof survive beyond fhe specfral gap closing. 

CONCLUSIONS 

The presenf work deal! exclusively wifh fhe bulk in- 
varianfs, while fhe invarianfs formulated in Ref. Ill can 
be regarded as boundary invarianfs since fhey are com- 
pufed from fhe boundary sfafes. We can already foresee 
a cormecfion befween fhese invarianfs, which we would 
like fo skefch briefly. The bulk-boundary principle de¬ 
veloped in Ref IISTl provides an equalify befween fhe 
bulk firsf Chem number and a cerfain specfral flow of fhe 
boundary sfafes. The laffer has a similar self-averaging 


properfy as fhe firsf Chem number. When applied fo 
fhe variafion of fhe elecfric polarizafion AP{y - Jy) of a 
1-dimensional sysfem wifh an edge, fhis bulk-boundary 
principle seems fo lead precisely fo fhe counfing of fhe 
specfral feafures of fhe edge sfafes performed in Ref. HJ. 
If will definitely be inferesfing fo make fhis connecfion 
more precise and see if in acfion for concrete models. 

The principle described in fhe presenf work seems 
fo apply fo any symmefry-sfabilized fopological invari¬ 
anf which can be formulafed in a real space represenfa- 
fion. Unforfunafely, fhere are imporfanf insfances where 
a real space represenfafion is nof yef available, such as 
fhe Kane-Mele Z 2 invarianf 1321 or fhe bulk fopological 
invarianfs for poinf-symmefry sfabilized fopological in¬ 
sulafors l33]| . For fhis reason, we have nofhing fo say 
abouf fhese invarianfs af fhis momenf, buf we hope our 
findings will spur a renewed efforf in fhis direcfion. Nev- 
erfheless, fhe sfrafegy does apply fo fhe Loring-Hasfings 
invarianfs B341I35I , or fo fhe spin-Chern numbers 13611371 . 
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